function [ s ] = evaluateGradientQuadraturePoint( p, q, qn )
%EVALUATEQUADRATUREPOINT Summary of this function goes here
%   Detailed explanation goes here

% r = p - q;
% 
% rn = norm(r);
% 
% rg = r / rn;
% 
% qnc = [qn(2); -qn(1)];
% 
% M = [rg(1)*rg(1), rg(1)*rg(2); rg(2)*rg(1), rg(2)*rg(2)];
% N = [0, dot(rg, qnc); -dot(rg, qnc), 0];
% I = diag(ones(2, 1));
% 
% prpn = dot(r, qn) / rn;

K = getKernelGradient(p, q);

s = tensorVectorDot(K, qn);
% s = tensorVectorContraction(K, qn);

end

